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The quark's behaviour while influenced by a strong stochastic gluon field is analyzed. An approxi- 
mate procedure for calculating the effective Hamiltonian is developed and the corresponding ground 
state within the Hartree-Fock-Bogolyubov approach is found. The comparative analysis of various 
Hamiltonian models is given and transition to the chiral limit in the Keldysh model is discussed in 
detail. 
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Nowadays we know that the mixing of the zero modes 
is the microscopic mechanism of the spontaneous break- 
down of chiral symmetry in the instanton hquid model 
[l[ . In this approach the quarks are considered in a given 
gluon background and the spectrum of the respective 
Dirac operator is calculated in order to be accompanied 
then by averaging over the gluon ensemble. It is believed 
that at low energy the zero modes are effectually over- 
lapped and the eigenvalues of Dirac operator spread over 
some range of virtualities. In other words, studying the 
behaviour of a single quark in external (stochastic) field 
one endeavours to guess the corresponding one-particle 
Green function but for the quark ensemble even now. 
Being unable to carry out this transition in detail from 
the first principles one is forced to devise the suitable 
approximations argued by some general theoretical rea- 
sonings. Also great care is needed in order to obtain 
proper thermodynamical limit with nonzero chiral con- 
densate. A lot of that happens to be in striking con- 
trast to the Nambu-Jona-Lasinio (NJL) model [2,] which 
is cognate to the instanton liquid model based actually 
on the similar multi-fermion interaction. Superficially, 
the main distinction consists in the appearance of some 
non-local formfactors instead of corresponding coupling 
constant. As to the microscopic consideration, the gen- 
eration of dynamical quark mass in the NJL model is 
caused by the reconstruction of the Hamiltonian ground 
state and the quarks manifest themselves already as the 
quasi-particles although the multi-fermion attractive 
force should be strong enough, roughly speaking. In this 
paper we emphasize that an instanton model and several 
other models which are based on treating the stochastic 
ensemble of strong gluon field become practically identi- 
cal in many aspects to the NJL model. 

Such an approach is motivated by the conceptual idea 
of an intricate nature of the QCD vacuum |4| having 
populated by intensive stochastic gluon fields of nontriv- 
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ial topological structure. Moreover, studying the corre- 
sponding cooled lattice configurations gives evidence of 
this component presence [5| and using the instantons in 
the singular gauge to fit the data turns out to be very 
fruitful 6] and allows one to evaluate the ensemble den- 
sity (around one topological charge per fm*) and the 
characteristic size of a saturating configuration (about 
1 GeV~^). Both estimates are in fairly good agreement 
with the corresponding results of instanton liquid model 
Q- Nevertheless, the keen search of various confining 
configurations is still going on Q , Q , [l3| in parallel with 
collecting the convincing evidences that the construction 
of self-consistent ensemble of such configurations is a too 
complicated problem (see, for example, the estimate for 
the (anti-)instanton ensemble done in Ref. [HI). 

Supposing the high-frequency component of stochastic 
ensemble of gluon fields as the dominating contribution, 
we develop, in fact, an effective theory[37| with apply- 
ing the procedure of simplified (averaged in time) system 
description which is widely used at studying the dynam- 
ical systems. Developing the effective theories which are 
discussed here has been launched to a considerable ex- 
tent by studying the behaviour of light quarks in the 
instanton gas (liquid) (l^ . The zero mode approxima- 
tion has provided for the quantitative picture of sponta- 
neous chiral symmetry breaking 13]. However, an effec- 
tive Lagrangian of the NJL type was soon received in Ref. 
fl^ by the direct summation of certain leading diagrams 
and the obtained vertices of multi-quark interactions oc- 
curred rather different from those calculated in the zero 
mode approximation. Analysis of heavy quark systems 
behaviour affected by the stochastic gluon fields has 
demonstrated that at constructing the respective effec- 
tive theory the cluster decomposition of generating func- 
tional [16] can be very efficient tool and the specific role 
of various characteristic correlation times has been clar- 
ified to classify the descriptions. These results together 
with the criticism of zero mode approximation [17] have 
contributed to widening the cluster decomposition ap- 
plications. This approach has been used to analyse the 
light quark behaviour plsj and it is interesting to note 
the effective Lagrangian has agreed with that obtained 
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in [l^. In the context of our interest here the cluster 
decomposition is called upon to describe the correlations 
in quantum system inspired by an external process. In 
this situation, as a matter of fact, the description of sys- 
tem behaviour is executed by averaging the generating 
functional. However, as we show in this paper such a 
procedure applied to the quantum system could be in- 
complete and it is more appropriate to base an analysis 
on the corresponding density matrix. Nonetheless, we 
argue here that in the 'white noise' limit (when the time 
intervals of stochastic impulses are very short) the proce- 
dure of averaging the generating functional occurs quite 
adequate. 

The form of the effective Hamiltonian obtained urges 
us to search the system ground state as the Bogolyubov 
trail function. The corresponding dressing transforma- 
tion will be analysed for various formfactors of effective 
Hamiltonian. In such an approach the quarks are al- 
ready treated as the quasi-particles and rather practi- 
cal way to get beyond the zero mode approximation ap- 
pears. It grounds on the method of simple iterations of 
corresponding integral equation solutions for the dress- 
ing transformation which quite stable unlike the results 
of mean field approximation , [13 ■ The different en- 
sembles are examined and their selection is stipulated by 
the requirement that one of their asymptotic forms would 
be the NJL model which plays a calibrating role in our 
calculations. The chiral limit of the Keldysh model with 
the correlator behaving as a (5-function in the momen- 
tum space is studied in detail and the singular behaviour 
of the corresponding mean energy functional is demon- 
strated. 



I. THE HARTREE-FOCK-BOGOLYUBOV 
APPROXIMATION 

We consider the quark (anti-quark) ensemble in the 
background of strong stochastic gluon field and suppose 

I 



here the summation over index s which stands to describe 
two quark spin polarizations and index c which should 
play the similar role for a colour is implied. Further we 
make concrete the form of the Dirac conjugated spinor. 
Fixing a spin polarization as it is known can be done by 
imposing an additional constraint on spinor (see, below). 
However, there is no direct analogy with the colour po- 
larization and the particular state should be fixed by the 
corresponding complete set of diagonal operators which 
includes the Casimir operators as well. In fact, this com- 



this field is so strong that we could neglect the gluon 
interchanging processes (quenched approximation). The 
stochastic gluon field is characterized by a correlation 
function and its particular form will be discussed and 
fixed below. The Lagrangian density is the following 

Ce = q iniiDfj, + im) q , (1) 

here q, q are the quark and anti-quark fields with covari- 
ant derivative = d^ — igA'^^t'^ where is the gluon 
field, t"' = \°- /2 are the generators of colour gauge group 
SU{Nc) and m is the current quark mass, /i — 1,2,3,4. 
We work in the context of the Euclidean field theory 
and 7^ mean the Hermitian Dirac matrices (7^ = 7^, 
{ifnlv} — 2 Sfj_^) in the chiral representation. Then the 
corresponding Hamiltonian description results from 



7i ^ nq - Ce , tt = = iq'^ , (2) 
and, in particular, for the nonintcracting fields we have 

Hq = —q {ij\7 + im) q . (3) 

In the Schrodinger representation the quark field evolu- 
tion is determined by the equation for the quark proba- 
bility amplitude 5" as 

* = -H^ , (4) 

and the creation and annihilation operators of quarks 
and anti-quarks a"*", a, b~^,b have no 'time' dependence 
and consequently look like 



(5) 

I 

plete definition of the spinor colour state is unnecessary 
for us here. All observables are usually expressed by sum- 
ming up the polarization states of some bilinear spinor 
combinations as the singlet and octet states and the sin- 
glet component is obviously playing the specific role. 

The density of interaction Hamiltonian can be pre- 
sented as 



Vs = 9(x) i"7^A;:(t,x) g(x) . (6) 



gQj(x) 



dp 



(2^)3 (2b4|)i/2 



[ a(p, s, c) Ua,(p, s, c) e^P" -t- 6+(p, s, c) v^^ip, s, c) e 'P"^] 
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The obvious dependence on 'time' in this Hamiltonian is 
present in the gluon field only. As it is mentioned above 
we are planning to work with the stochastic gluon field 
implying the random process for which one may define 
only a probability of realizing some gluon configuration. 
Such a nature of gluon field urges (and allows) us to de- 
velop the approximate procedure for describing the quark 
field treating (|3]) as a probabilistic process. Then the sys- 
tem states are described by the corresponding averages 
(over a 'time' or an ensemble according to the ergodic hy- 
pothesis). However, in the quantum theory we face one 
difficulty in this way because ^' is a probability amplitude 
and an immediate averaging of (^) can be insignificant. 

Studying a mean probability density (\I>vI/) looks more 
promising and can be realized by complicating the pro- 
cedure of continual integration |19j . In Appendix I we 
analyse convincing quantum mechanical example to il- 
lustrate the difference between two approaches. One of 
those is based on constructing the corresponding density 

matrix (vj/vj/), and the second approach does use the rele- 
vant averaging of the functional ('5). We argue the latter 
could be practical for application in the 'white noise' limit 
with the (5-like time correlation function. Adapting these 
ideas to the gauge theories we should obviously strive to 
operate with the gauge invariant quantities which include 
an ordered exponential, at least. Unfortunately, such a 
program in what concerns the ensemble consideration is 
still very far to be realized. However, it is clear that ap- 
plying the averaging procedure would result in putting 
in an appearance of a set of corresponding correlation 
functions {A^), {A'^) etc.jsl 

In the interaction representation, where ^' = e-^o*$, 
Eq.(|4]) can be rewritten as 

$ = -y$ , V = e^oVse--^"* . (7) 

Now the 'time' dependence appears in quark operators 
as well. Now we remind some features of the averaging 
description as formulated in Ref. Presenting Eq.([7]) 
in the integral form as 

$(i) = $(0) - I dr V{t) $(r) , (8) 
Jo 

where $(0) is an arbitrary initial state of ensemble and 
performing another iteration one receives 

$(t) = $(0)- [ dTF(T)$(0)+ / dTV{T)[ dT'y(r')$(r'). 
Jo Jo Jo 

By averaging the fast-changing component and uncou- 
pling the correlators one approximately approaches the 
long wavelength component ($) in the highest order (also 
taking into account that {V) — 0) as follows 

($(t)) « $(0) + fdT /V(V^(t)V^(t'))($(t')). (9) 
Jo Jo 

Certainly, it is assumed the characteristic correlation 
time of stochastic process is smaller than the time char- 
acteristic for the process ($). By differentiating Eq.® 



it is easy to get rid of the initial condition $(0) and to 
have 

im) =+ f dr' {V{t)V{T')) (<i>(r')) . 
Jo 

Actually this equation should describe a steady-state pro- 
cess and at reversing a time the solution, in general case, 
will not return to the initial magnitude $(0). Changing 
the integration variable as r' = t — t one comes to 

($)=+/" dr {V{t)V{t-T)) {<i>{t-T)) . (10) 
Jo 

It is usually supposed the correlations are quickly decay- 
ing then the upper limit of integration might be changed 
for oo and in order to deal with the local process it is 
well justified (without a precision loss) to change the ar- 
gument of function ($) for t. Eventually, as a result we 
have 

/>CXD 

m))=+ dr{V{t)V{t-r)){m)- (11) 
Jo 

(The requirements to validate the factorization of the 
long-wavelength component are discussed, for example, 
in [16].) Implementing the approximation (jlip in the 
quantum field theory models, we run into the trouble at 
trying to get the most general form of correlation function 
if the characteristic quark and gluon correlation times are 
comparable. Fortunately, if the quark fields are consid- 
ered to be practically constant on the gluon background 
the problem receives essential simplification. The gluon 
field contribution may be factorized as a corresponding 
correlation function {A'^^{x)A^^{y)) [23|. Recent lattice 
measurements of this correlation function provide us with 
a reasonable arguments to interpret the result as gluon 
'mass' generation 300 - 400 MeV) in the momentum 
region of order 200 MeV [IH. 

It is curious to notice that the averaging over ensemble 
('time') in the right hand side of Eg. pTj) is performed in 
both the correlator and ($(i)). It means that by resum- 
ming and averaging a certain class of diagrams in the 
quantum field theory models, one may take into account 
high order correlator contributions in different ways if 
the form of function ($(i)) is specified. Besides, the cor- 
relation functions in models interesting to us should be 
translation invariant and it implies that the correlator in 
Eq.pT]) has the following form 

{V{t)V{t-T))=F{T) , 

i.e., for example, an one-dimensional process after hav- 
ing done the integration in Eq. (|lip will be described by a 
constant which characterizes the slow process. In quan- 
tum field theory for the problem we are interested in, the 
correlator connecting two space points 
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dx g(x, t) g(x, i) dr / dy g(y, i ~ r) ^"7, q(y, t ~ r) ^'(^^(i, x)A^(i - r, y)) ($(<)) 



r 



appears instead of a constant. Assuming the correlation form the inverse transformation to the Schrodinger repre- 
function is rapidly decreasing with time we change the sentation. Then introducing the function x — e~^''*($) 
'time' t — T dependence in the quark fields for t and per- we have[^] the following equation 



X = -H^nd X , 



im)q- q ^7,, q j dy q' t'j, q' dr g^A^A^) 



(12) 



r 



where q = g(x), q = q{x), q' = q{y), q' = q{y), = 

A;(t,x) and A',b ^Al{t~T,y). 

In order to receive the final result we should fix the 
form of correlation function. In this paper we rely on the 
stochastic ensemble of (anti-)instantons in the singular 
gauge 



N 

E 

1=1 



(13) 



and instanton solution reads as 



If 2 

<P{q) = — K2{qp) - 

q^ \ q^p^ 



(27r)4 ^^"^^ 



(14) 



where is the modified Besscl function of imaginary 
argument, p is the instanton size, the matrix lo appoints 
the pseudo-particle orientation in colour space, z is the 
coordinate of instanton center and f\ stands for the 't 
Hooft symbol. The distribution of the pseudo-particle 
orientation in colour space is supposed to be homoge- 
neous ^ duj as well as the probability to find a pseudo- 
particle in the volume element is proportional ~ dz/V 
where V is the volume of the system under consider- 
ation. Apparently, specifying the saturating configura- 
tion in the form of Eg. p^ is, in a direct way, the gauge 
fixing procedure. Calculating the quantum corrections 
for every single pseudo-particle in one-loop approxima- 
tion (what corresponds to the zeroth order of the N/V- 
expansion), and exploiting the variation principle 0], [lH 
allows one to ascertain the size distribution of pseudo- 
particles. In this way it is possible to attach clear mean- 
ing to the functional and to construct in the thermody- 
namical limit limy^oo N /V — > n the state possessing a 
negative energy density and developing a non-zero gluon 
condensate. (Uncertain interrelation of perturbative and 



non-perturbative contributions into the path integral [22 | 
makes the computability of generating functional highly 
nontrivial as for now.) In Eq. (|12p we imply the corre- 
lation function integrated over the 'time' for which we 
receive in the highest order in the density n of (anti- 
)instanton ensemble 



dx4 {Al{x)Al{y)) 



dx, {Al{x)Al{y)) = 



4(47r^)^ Sab np^ . . w 

^ e^P(^-«) 0(-p)0(p) \ 2^ 5{p^) . 

The first equality is valid due to the symmetry properties 
of instanton solution. Then the correlation function can 
be presented as 



(47r2)2np4 2 5ab 



[I{p) 5i_iu - J^i,^{p)] 



52 - 1 

lip) = H''P)<Pip) > Jtjip) = PiPj 0(-p)0(p) 



(15) 



We suppose in what follows the various stochastic en- 
sembles of gluon fields are characterized by their pro- 
file functions /(p), J^v{p) and analyze the contribution 
of quadratic correlator only. However, this deficiency 
of fixing the gauge implicitly for the truncated system 
is compensated, in a sense, by our investigation of full 
spectrum of reasonable correlation functions (including 
an opposite limiting correlators when they are extrap- 
olated even into the perturbative region). Recent con- 
siderable progress in studying the confining configura- 
tions of lattice gauge theories, in particular, revealing 
the monopole clusters and their role in confinement (see, 
for review [23[) as well as detecting the specific features 
of quark behaviour in the uncooled configurations and 
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the indications that low-lying Dirac eigenmodes are lo- 
calized on the objects of dimension inherent in monopoles 
and vortices [l^ looks entirely urging. But these results 
are also bringing the perilous tendencies because reveal 
the some features of lattice gauge theories common with 
the compact electrodynamics. This fact makes almost 
inevitable to draw in the singular (in the continual limit) 
objects and to give them an underlying physical meaning. 
Searching the formfactors (the corresponding ensembles 
of saturating configurations) interesting for applications 
one should compare to the reasonable results for the four 
nonets of light mesons obtained in the NJL model. Ap- 
parently, the constants of effective four-quark Hamilto- 
nian (scalar, pseudo-scalar, vector and axial-vector chan- 



nels) and the parameters of integral saturation (cut-off) 
should be comparable. Seems, it might be carefully sup- 
posed that those singular objects (at still an unknown 
scale) have to reproduce the major features of successful 
effective NJL Hamiltonian after the corresponding aver- 
aging. In any case, the problems to find the specific fea- 
tures of such singular objects which admit their experi- 
mental identification and to analyse the quark behaviour 
in the ensembles of monopoles or vortices are of really 
great interest [25l |. 

With such a form of induced four-fermion inter- 
action we are going to search the ground state as 
the Bogolyubov probe function with vacuum quantum 
numbers [Iq [H, ^ 



\a) ^ T |0) , 

T = Tlp^s,c exp{ ip [ a+(p, s, c) 6+(-p, s, c) -I- a(p, s, c) &(-p, s, c) ] } 

I 



(16) 



which is defined by minimizing mean energy 
E=(a\H\a) , 



(17) 



tion operators of quasi-particles {T i = for fcrmions) 
A = TaT-^, B+ = Tb+T-^, 



here ip = f(p) and |0) is the vacuum of free Hamiltonian, 

i.e. a(p, s, c) |0) = 0, b{p, s, c) |0) = 0. Introducing with we present the operator Eq.([5|) as, with the Dirac conju- 
the dressing T transformation the creation and annihila- gate spinor 

I 



9(x) 



dp 



(2^)3 {2\p,\y/^ 

dp 1 

(2^)3 {2\p,\y/^ 



[ A{p, s, c) U{p, s, c) e'P'' + B+{p, s, c) V{p, s, c) e~'P^] 
[ .4+(p, s, c) C7(p, s, c) e-^^- + i3(p, 5, c) F(p, s, c) e^P-] 



(18) 



where the spinors U and V are defined as 

U{p,s,c) — cos((^) u(p, s, c) — sin(c/j) «(— p, s, c) 
y(p,s,c) = sin((^) u(-p, s, c) + cos(i^) u(p, s, c) 



(19) 



with U{p,s,c) = C/+(p, s,c) 74 and V{p,s,c) = 
l^+(p, s,c) 74. Now we have to specify the choice of 
spinors in the Euclidean variables. They obey the Dirac 
equations 

(j) — im) u{p, s) — , {p + im) v{p, s) = , (20) 

(with p — piji + P7) and additional constraint which 
fixes the spinor polarization 

175 s u{p, s) = u{p, s) , 175 s v{p, s) = v{p, s) , (21) 

where 75 = —71727374, and the four- vector s is nor- 
malized to unit and orthogonal to the four-vector p, i.e. 



1, (ps) = 0. It could be, for example. 



Si 



(pn) 



(pn) p 



im (p4 — im) 



where n is an arbitrary unit vector. If the covariant nor- 
malization conditions are satisfied 



uu — 2im 



vv — —2im , 



(22) 



the spinors are defined with the precision up to phase 
factor. All these conditions allow us to formulate the 
following matrix representation 



P + iiTi 

u[p, s) u{p, s) = — - — (1 + 175 s), 
v{p, s) v{p, s) = — - — (1 -I- 175 s) . 



(23) 
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Calculating the mean energy Eq. (fT7|l we meet spinors 
with opposite moments. We introduce the four-vector 
q = {pi, — p) in order to simplify notations. Using the 
projection operator we can express the spinor v{q, s) 
through the spinor u{p,s) (see '28]) 

v[q,s)^a—-. u{p,s) . (24 

—2im 2 



The coefficient a is fixed by the covariant normalization 
Eg. ip^ up to the phase factor as 

* 2 rn? m 

aa— -- — 7T — , a = 1 — r ■ 

[pq) + |p| 

Then summing up over the spinor states results in 



Vu(g,s) i;(p,s) = a ^ + (p - im) , V w(p, s) u(g, s) = a (p-im) 

^ — ^ /im ^ — ^ 



2im 



q + im 
2im ' 



''^^u{p, s) v{q, s) — a {p + im) 



q — im 
2im 



^v{q,s) u{p,s) 



q — im ^ 

a (jp + im) 

2im 



(25) 



The polarization in which the momentum p and unit po- 
larization vector n are orthogonal (pn) = turns out to 
be the most convenient for handling. In such a situation 
both operators p and q commute with 75 s and the polar- 
ization directions of quark and anti-quark could be taken 
identical (although in general case they should be two 
different directions). Then the summation over polariza- 
tion of quarks and anti-quarks is performed separately in 
the final equations. It allows us not to control the oblig- 
atory constraint to have the vacuum quantum numbers 
of the pairs present in the intermediate calculations. 

When calculating the mean energy Ea. (|17[) nontrivial 
contribution bilinear in quark operators comes from the 



terms of type {cr\B B^\a) (remember B\cr) = 0, A\a) = 
0). The contribution from the terms like {<j\A A^\a) is 
absent because of the particular representation of bilocal 
operator we are using as qq (then quadratic terms are ex- 
pressed by the spinors V, V . Due to the similar reasons 
the four-quark operators develop only two nonzero con- 
tributions {a\B B+ B' B'+\a) and {a\B A a' + B'+\a). 
The first combination corresponds to the contribution of 
so-called tadpole diagrams and the latter is related to 
the asterisk[4l| diagrams. As a result the four-fermion 
interaction contribution can be presented in the follow- 
ing form 



+ J ^IpIwp'^i F„.(p,5,c) tf^^^pVsi{p,s,c) U,k{p\s',c') tin';5Upj{p\s',c'){A^^{p + p')). 



Here (^^A'''^) is the Fourier transform of the gluon cor- 
relator and the summation over spinor and colour in- 
dices is implied. The contribution of the first tadpole 
diagram is an identical zero due to completeness of the 
spinor basis in color space, giving a unit color matrix 
(in particular it is valid for colour singlet quark config- 
uration). In electrodynamics the considered term pro- 
vides a dominant contribution. But it is interesting to re- 
mark that the singular character of photon propagator in 
the infrared region makes this abclian theory even more 



complicated to research than in the nonabelian one. In 
the compact U{1) electrodynamics (on a lattice) the in- 
frared behaviour of correlation function is formed by the 
monopole contributions but nowadays it is still difficult 
to define a scale where these effects show up themselves. 
In the octet channel of nonabelian theory we obtain the 
quark repulsion ~ — l/(47Vc) and therefore this regime 
might be omitted when searching the minimum of mean 
energy Eq.([T7|). Then for the spinors with polarizations 
summed up we have 
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VV = P474 + cos{9) (p7 — im) 



a + a / X / 9 

— sm[9) (p — im pj) 

2im 



— a+a /x/9 

= P474 + cos{9) (p7 + im) + /^-^ 



2im 



sm{6) {p^ + im P7) , 



where angle 9 = 2(p. In the formulae above the phase 
inherent in the sum a + a (a spinor is defined up to 
such a phase) is still indefinite. The direct analysis of 
the mean energy functional demonstrates that the most 
preferable value of the phase factor (responsible for the 
colour interaction of quarks) is the value when the coef- 
ficient a appears to be a real number. For definiteness 



we put a = +\m\/p. The curious fact is that the results 
of summation are not equal {VV{m) — UU{—m)) and 
they coincide in the chiral limit m — Q only, i.e. particles 
and antiparticles formally generate the different contri- 
butions. 

The direct calculations lead to the following result for 
the mean energy (flT)) 



{(T\Hind\'j) = - 

dpdq ( 



-G 



dp 2Nc pI 

P4 94 



(2^)3 



(27r)6 \" 



{31- J) 



(2/4 



b4||'74| 

b4||g4| 



(1- cos 61) 
+ (4/ - J) 



P g 

b4||g4| 



cos ( 



p 

m 



m 

— cos ( 

q 



■ sm ( 



(26) 



here we designated p = |p|, q = |q|, / = I{p + q), = 

Jtj{p + q)^, = ELi -^"^ pI+p^ = ql + = 

6' = 9{q), G = (47r^)^ np^ and as a matter of convenience 
we singled out the colour factor G' = j^'i_^ G. To obtain 
this result we performed the regularization (subtracting 
the free Hamiltonian Hq). It results in the presence of a 
unit (together with — cos 9) in the parentheses of the first 
integral. Let us also remind that in the Euclidean space 
p\ is a negative magnitude. Then with Eq.([26|) available 
one can find the most advantageous value of the angle 9 
from the condition 

'^^^\H^na\cT) ^ Q ^27) 
d9 

Henceforth we characterize the different stochastic en- 
sembles of the gluon fields by their profile functions I{p), 



J^u[p)- 



II. NAMBU-JONA-LASINIO MODEL 

Now let us consider the example in which the corre- 
lation function behaves in the coordinate space as the 
(5-function (simply we assume J^u{p) — 0)- Actually, it 
corresponds to the Nambu- Jona-Lasinio (NJL) model [2] ■ 
As well known the regularization is required to obtain an 
intelligent result in this model. We adjust the NJL model 
with the parameter set given by Ref. '29*1, and limit the 
integration interval over momentum in Eq. (|26p with the 
quantity |p| < A (A = 631 MeV). Then the functional 
(j26p is written in the following form (unessential terms 
contributing the constant values are omitted) 



W ^ 



dp 

(2^ 



b4|(i 



COS( 



Gt 



|P4| 



m 

sin 9 cos ( 

^ P 



(2^)3 \q,\ 



' . m 
sin 9 cos ( 



(28) 



here m = 5.5 MeV. The equation to calculate the optimal where 

angle 9 (l27l) reads as ^ 

(2^b? V ^"7^°^ / ■ ^^^^ 
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The constant of four-fermion interaction is G = -^^G 
while expressed in the same units as the mean energy 
functional in Eq.(l26|). For the NJL model Eq.([29l) makes 
it possible to contract a functional space in which the 
minimum of mean energy functional can be realized. This 
equation parameterizes the function d{p) on the whole 
interval p £ [0, A] of searching the solution. Moreover, 
Eq. ([M)) itself does not impose any restrictions on the pa- 
rameter M which may be any real number. Then the 
functional (|28p simply becomes the function of param- 
eter W{M). Now if one expresses the trigonometrical 
functions via parameter M it is possible to make the 
representation of minimizing function and the result of 
its integration (130]) agree. As a result we receive three ex- 
tremal points, two of them correspond the minimal points 
with negative and positive values of A'l and the negative 
value conforms to the state of more stability. The point 
of unstable equilibrium is located in the vicinity of co- 
ordinate origin ~ m. The induced quark mass for the 
parameter magnitudes fixed is M = —335 MeV and the 
quark condensate 

(cr|gq|cr) = — — / dp - — - (p sin 6* - m cos 6*) , (31) 

J a \P4\ 

develops the magnitude of {(j\qq\(j) = —i (245 MeV)^. 
The characteristic constant of the four-fermion interac- 
tion is equal to G'/(27r^) = 1.34. In what follows we rely 
on these quantitative results. 

The situation, if the correlator J^v{p) is not equal to 
zero and has the same form of the (5-function in coor- 



dinate space, can be similarly analyzed. The numeri- 
cal analysis done teaches that its influence can be es- 
sential but we do not show these results due to the 
lack of any phenomenological estimates of the correla- 
tion function magnitude. The non-local version of the 
NJL model in which the correlator has the separable form 
/(p, q) = K{p) K{q) can be similarly analyzed. In fact, 
it again displays the above mentioned property which re- 
places the functional analysis for the analysis of function 
dependence on some parameter. Although one important 
difference does exist and it shows that the procedure of 
integral cutting off is unnecessary for the functions K{p). 
The regularization is naturally performed by the K{p) 
kernel and so strong regularization is caused by the sepa- 
rable form interaction kernel. Certainly, such a property 
can manifest itself in much weaker form for more realistic 
correlators. 



III. THE KELDYSH MODEL 

Here we are going to analyse the limit in which the cor- 
relation function has a 5-function form in the momentum 
space 

/(p) = (27r)3 G 5(p) . 

This limit is an analogue of the Keldysh model which is 
well known in the physics of condensed matter 'so'l and 
the mean energy functional (|26[) develops the following 
form in this case,42j 



W{m) 



dp 

(2^ 



\pi\ (l-cos6')-G 



b4P 



sin 9 



■ cos 9 



(32) 



The optimal values of angle 9 arc determined by the so- 
lutions of the following equation 

|p4p sin6'— 2G (pcos6' -f- msin^) (psin6' — mcos0) = 

(33) 

and we start analyzing these solutions in the chiral limit 
m = 0. One of the solutions corresponds to the zero 
angle 9 = but the non-trivial one takes the form 

cos0=^. (34) 

Both the positive and negative angles 9 are suitable as 
the solutions because of the parity (positive) property of 
the functional ([5^ and these (real) solutions (additional 
to the trivial one) exist on the limited momentum inter- 
val p < 2G. There are one real solution for the trivial 
angle and two imaginary (complex-conjugate) solutions 
beyond this interval. Analyzing the NJL model above 



we noticed its very convenient property when the solu- 
tion 9{p) is defined on the whole interval and, in fact, 
the functional is parametrized by a single number which 
is the integral M. In the Keldysh model the situation 
is much more sophisticated and the phase portrait of 
its solutions in the chiral limit (for example, sin0 as a 
function of momentum p) consists of two arches (with 
positive and negative sm9, see Fig.l) and straight line 
corresponding to the trivial solution. Thus, the semi- 
axis p € [0,oo) can be divided into two parts. There are 
three branches (solutions) at the interval p £ [0, 2G], two 
of those correspond to the positive and negative angles 
9 ^ ± arccos(p/2G) and the trivial one where 9 = 0. At 
the interval p £ [2G, oo) only one trivial solution 9 = 
exists and in order to construct the solution on the whole 
semi-axis p £ [0, oo) one has to add the trivial solution 
on the interval p £ [2G, oo) to any detached branch of 
solutions on the interval p £ [0, 2G] . It is easy to see that 
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we obtain for the mass Mg which characterize the angle 
at the top arch the following result 



100 200 300 400 500 

p MeV 



FIG. 1: Phase portrait of the Keldysh model, sinO as a func- 
tion of momentum p(MeV). The dotted curves correspond to 
the solutions in the chiral limit m — 0. 



making use the imaginary branches of solutions leads to 
the significant growth of energy and just because of this 
fact they are uninteresting. The other potentially inter- 
esting functions 9{p) for which it is reasonable to search 
the functional minima could be received if the interval 
p e [0, 2G] is subdivided into smaller intervals and then 
for each interval when continuing the function to the next 
interval (for example, to the direction of the momentum 
p increasing) to use two other branches as well as of the 
results of continuation on the same branch. Apparently, 
it results in the piecewise continuous function and un- 
like the NJL model here we have no parameter which 
restricts the function and watches its integral character- 
istics. In the chiral limit all the solutions (trajectories) 
constructed in such a way will acquire strictly fixed (fi- 
nite) value of the functional W{0) (it will be observed 
that the functional docs not contain the derivatives of 
angle in momentum). For example, the trajectory which 
is going along the top arch at the interval p G [0, 2G] 
and continuing longer as a trivial solution to the whole 
semi-axis leads to the magnitude 

157r^ 

(similarly for the top negative arch). The chiral conden- 
sate (ISTI) turns out then to be 



{a\qq\a){0) = 



2^ 



(and for the solution along the negative arch we have the 
opposite sign). The mean energy and chiral condensate 
equal to zero for the trivial solution, i.e. (Wo(0) = 0, 
{(j\qq\a)o{0) = 0). Clearly, these piecewise continuous 
functions will lead to the magnitudes of functional W{0) 
which fill up the interval [Wo(0), VF±(0)] densely, a sim- 
ilar pattern takes place for the chiral condensate. With 
the natural parametrization 



Mg = (4G^ 



P 



.1/2 



(36) 



Mg 



smt 



(p2 + Mgy/^ 



(35) 



It is interesting to notice that then the respective en- 
ergy of nontrivial solutions E{p) = ^Jp"^ + Mg becomes 
constant E{p) — 2G. 

After having done the analysis in the chiral limit which 
is shown by the dotted lines in Fig. [l]we would like to 
comment on the situation beyond this limit, i.e. where 
m ^ 0. The evolution of corresponding branches is avail- 
able on the same plot[T]where the behaviour of 9{p) as the 
function of momentum p in MeV is shown for the solution 
of Eq. (j33p . The semi-axis p G [0, oo) where we are search- 
ing the solution can be subdivided into two sectors which 
are demonstrated by the vertical dashed line on the plot. 
Three solutions denoted by a, 6 and c are developing at 
the first sector denoted in Fig. [1] by I. Besides, there 
are three solutions at the second sector denoted by II, 
one real solution designated as A for the negative pairing 
angle and two complex-conjugate roots with the positive 
real parts. The imaginary parts of solutions are plotted in 
Fig. [T]by the dashed lines. The solution A in the domain 
II develops the behaviour of 6* ~ —^^p- with increasing 
momentum. As in the chiral limit the minimum of the 
mean energy functional W{m) can be obtained with the 
piecewise continuous functions which are properly repre- 
sented by the trajectories a A, bA, cA (for real solutions). 
The first symbol of this complicated designation implies 
the branches a, b, or c at the first sector, the second sym- 
bol corresponds to the branch at the sector II. Thus, at 
low momenta we start with the solution of branches a, b 
or c, then relevant solution passes to the branches inter- 
changing its position in any subinterval. But in any case 
there is only one way to continue the real solution when 
momenta goes to infinity and it is related with the branch 
A when the angle is going to the zero value. Moreover, 
if the angle 9 could take strictly zero value in the sector 
II then the second term of Eq. ([5^ leads to the singular 

2 o 

contribution coming from the term ^ cos 9 with the 
linear divergency at large momentum. Besides, the other 
terms develop the logarithmic divergencies as well. It is 
an amusing fact the mean energy functional out of chiral 
limit goes to an infinity at any nonzero value of current 
quark mass m although in the chiral limit VK(0) is well 
defined. (It is worthwhile to remember here the current 
mass singularity of zero mode approximation which was 
discovered in Ref. flT^). The same conclusion is valid for 
the chiral condensate (see Ea. ((3T|) in which the first and 
second terms are developing the linear and quadratic di- 
vergencies, respectively. We could conclude here that if 
the cut-off factor is not used in the integrals when deal- 
ing with the solutions on the whole axis the functional 
W{m) and quark condensate {a\qq\a){m) are ill-defined. 

Let us remember here that by definition the approxi- 
mation (jl2p should describe the quark behaviour in the 
background of stochastic gluon field (which is averaged) 
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FIG. 2: The equilibrium angle 9 (in degrees) as a function 
of momentum p in MeV. The solid line shows the result of 
the NJL model and the dashed line corresponds to the most 
stable branch of the Keldysh model, the current quark mass 
is taken as m = 5.5 MeV. 



at low energies. Then it looks quite natural to introduce 
an effective cut-off (in momentum) parameter A. The 
condition for factorization of gluon and quark field con- 
tributions gets broken at the momenta above A. In such 
a situation the dependence of mean energy and quark 
condensate on the current quark mass is defined not only 
by the form of integrand but by the value of parameter A 
as well. And if this value is pretty large A ^ A/ (where 
M is the dynamical quark mass) the dependence on m of 
all the observables is mainly defined by the magnitude of 
cut-off parameter A because of the singular character of 
integrals (for example, for the NJL model this magnitude 
could be estimated as ~ IGeV). Obviously, it means in 
order to get the dependence of observables on the current 
quark mass we need to draw essential additional informa- 
tion. 

As to the possible interpretation of the singularities 
available in the mean energy functional we could assume, 
for example, the mechanism similar to the Cooper pair- 
ing which takes place at every scale of the increasing 
momenta A, Ai, .... Certainly, we should correct the 
existing results about four-quark interaction potential to 
put the pairing effect on realistic ground. For example, 
the contribution of the stochastic configurations like the 
small size instantons which is exponentially suppressed 
is hardly relevant to provide an efficient pairing mech- 
anism for the momenta above A. Apparently, the hard 
gluon exchange looks like a more adequate mechanism 
at small distances. Then the gluon correlation function 
in Eq. (jl2p should be transformed in the corresponding 
gluon propagator. The effective four-quark interaction 
we are interested in can be derived by the quasi-average 
formalism [26| which approximates smoothly the infrared 
and ultraviolet momentum regions although an alterna- 
tive scenario could also be q uite meaningful (see, for in- 
stance, the discussion in [31[). The fact that the Cooper 
attraction is still large enough despite the coupling con- 



stant weakening could signal the dominance of more fun- 
damental fields at very small distances. (Here it is easy 
to see all the different models might be classified by the 
convergence of the integral over momentum with the con- 
stant of four-fermion interaction 

Ig = j dp G{p) , 

as the integrand. The model falls under the category of 
a singular one if this integral diverges.) 

In Fig. [5] we compare the equilibrium angles 9 for the 
NJL model (solid line) and the aA solution of the Keldysh 
model (dashed curve) as the functions of momentum p 
in MeV with the current quark mass m = 5.5 MeV. It 
is interesting to notice that out of the chiral limit the 
solution (which has a spherical symmetry) passed over 
zero at p = (see. Eg. ([57)1 ). Besides, Fig 12] demonstrates 
that out of the chiral limit the pairing process becomes 
essential not at zero momentum value (as it takes place 
in the chiral limit) but it is shifted to the magnitude 
about pg ~ 40 MeV for the fitting parameters used. For 
example in the NJL model it can be obtained 

pg — [m \M — mj]"'"^^ . 

The quantity rg = 1 /pg determines the characteristic size 
of the region which is efficient for the pairing process. 
In the chiral limit this region is formally extending to 
infinity. The curves shown in the plot correspond to the 
opposite, in a sense, limiting regimes and it is interesting 
to evaluate where the model with more realistic correlator 
could be found out. 

One of the important motivations to study the Keldysh 
model was the question of a natural regularization which 
presents for the interaction with separable kernel. We 
have seen that in the chiral limit for the kernel with most 
extensively expressed regularizing property as, for exam- 
ple, the momentum (5-function, both the mean energy 
and chiral condensate are well defined. Out of the chiral 
limit the unexpected singularity appears. In Ref. [32] it 
was discussed the possibility of continuing the mean en- 
ergy functional and the quark condensate by performing 
the respective regularization. As well known the me- 
son masses in the NJL model can be presented by the 
quark condensate what hints the corresponding expres- 
sions in the Keldysh model could be singular as well and 
one needs to perform another regularization to provide 
them with clear physical meaning. However, despite the 
present singularity of chiral condensate the meson ob- 
servables are finite and are well matched with the ex- 
perimental mass scale (see |3^). The reason to have 
these meson observables as the smooth functions of cur- 
rent quark mass is in the regularizing role of additional 
vertex formfactors which enter the meson mass formu- 
lae. Then we may summarize that it does not make sense 
to debate about an absolute value of quark condensate 
(in vacuum) for considered mechanism of spontaneous 
chiral symmetry breakdown because its magnitude de- 
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FIG. 3; The parameter Me (sohd and dashed Unas) and quark 
condensate (without an imaginary unit, in power 1/3) in MeV 
(shown by points calculated in the chiral limit m = 0) as a 
function of the constant G for the Gaussian correlator. The 
solid line is calculated with the current quark mass m = 5.5 
MeV and the dashed line is calculated in the chiral limit. 



pends on the particular observable (characteristic mo- 
mentum of saturation) which is used for extracting this 
data. In the Appendix II we compare the results ob- 
tained in the Hartree-Fock-Bogolyubov approach for the 
NJL and Keldysh models with the results of mean field 
approximation. 



IV. THE EXPONENTIAL AND GAUSSIAN 
CORRELATORS 



Here we turn to a more realistic situation and an- 
alyze the solutions possessing a spherical symmetry in 
the regime where the correlation function /(x) is rather 
quickly decreasing with the distance increasing. Perform- 
ing the integration over the azimuthal angles we can get 
the equation to derive the optimal angle in the following 
form 



J 



4G 

TT 



cos d ^ sm 

P 



dqdx - — r {qsvaO' — mcosO') I{x) sm{px) sm{qx) — 0. 



(37) 



Considering the solution behaviour at high momentum 
values p we are interested in analyzing solutions in which 
the angle 9 is going to zero magnitude. Assuming the 9 
value to be rather small we expand p7p up to the terms 
of the 9 order and have 



p^9 



AG r°° 

— dq{q9' 

TT Jo 



m) I{p, g) = 



If the function 9 is decreasing faster than 1 /q the most 
essential contribution to the integral comes from the term 
proportional to m and if the kernel J(p, q) is integrable 
the asymptotic behaviour has the following form 



AG m 

IT p"^ 



dq I{p, q) 



Let us consider now two concrete examples, with expo- 
nential behaviour of the correlator /(x) = exp(— a |x|), 
and with the Gaussian behaviour /(x) = exp(— x^). 
The integration over x can be performed exactly for both 
cases and the kernels of integral equations look like 



POO 

I dx e~°'^ sin(px) sm{qx) 
Jo 



1 



1 



2 \a'^ + {p-q)^ a'^ + {p + qy 



for the exponential correlator and as 

/ dx e~°' ^ sin(pa;) sin((7a;) — 
Jo 



Aa 



for the Gaussian one. Now let us hold the contribution of 
the first term only at large momentum values p for both 
examples. Then as a result the corresponding asymptotic 
behaviours are expressed by the constants which are de- 
fined by the integrals with the kernels I{p,q). It allows 
us to conclude that we have again the singular functional 
for the mean energy out of the chiral limit. The pa- 
rameter Mg (see Ea. ([35t and the quark condensate as 
functions of the constant G for the Gaussian correlator 
(both obtained by the numerical computation of Eg. ([57]) ) 
are depicted in Fig. O The solid line demonstrates the 
solution with the current quark mass m — 5.5 MeV and 
the dashed line is calculated in the chiral limit as the 
quark condensate presented by the points. The intrinsic 
change of the parameter Mg generation out of the chiral 
limit is easily seen. The similar features are observed for 
the exponential correlator as well. 

Unfortunately, it is a very serious problem to get all 
the solutions of the nonlinear integral equation ([57|) and 
here we are working with only one of its (the most sta- 
ble) branches. As it was demonstrated above such a sit- 
uation generates a lot of difficulties for extracting a re- 
liable information on the observables out of the chiral 



12 




^^.0 6.5 i .0 1.5 2.0 
P* P 



FIG. 4: The parameter Me in MeV as a function of momen- 
tum p in GeV which corresponds to the best fit of 'experi- 
mental data'. The solid curve is calculated for the Gaussian 
correlator with m = 5.5 MeV, the dashed line is calculated 
for the same correlator but in the chiral limit. 



limit. Due to this reason we calculate here the dynami- 
cal quark mass (M = —335 MeV) and chiral condensate 
{\{a\qq\a) \ = (245 MeV)^) in the chiral limit collating the 
dynamical quark mass with Mg and fitting the parame- 
ters a and G. The parameter a for the exponential and 
Gaussian correlators reads as 

fle^ = 0.15 GeV, ags = 0.16 GeV . 

The most suitable values of G are equal to 

Ge. = 0.35 , Ggs = 0.31 

and |M|^| = 338MeV, \{a\qq\(j)f' \ = (228 MeV)^, 
|Mf I = 340MeV, |((T|g(7|cr)f | = (245 MeWf. Actually 
we can collate the gotten value of the four-fermion inter- 
action constant with the packing fraction parameter of 
which the basic one for the instanton vacuum model is 
G = (47r^)^ np^. The result of this exercise np^ ~ lO^'' 
is quite realistic. However, there is a pretty serious dis- 
crepancy in the estimates of characteristic configuration 
size (we should keep in mind the calculations are done in 
the chiral limit). 

The parameter Mq as a function of momentum p in 
GeV calculated with the parameters corresponding the 
best fit of 'experimental data' is depicted in Fig. 21 The 
solid curve is obtained for the Gaussian correlator with 
TO = 5.5 MeV and the dashed line is calculated for the 
same correlator but in the chiral limit. Here we do not 
mention the results obtained for the behaving exponen- 
tially correlator because they practically coincide with 
the results for the Gaussian correlator. The parameter 
pe is estimated at the current quark mass to = 5.5 MeV 
to be as pe ~ 150 MeV, i.e. of the 7r-meson mass order. 
The treatment of the correlator with instanton profile 
together with the detailed analysis of exponential and 
Gaussian correlators is worthy of special paper and will 
be studied in the next paper. 



V. CONCLUSION 

In the present paper we undertake the efforts to sys- 
tematically study the quark's behaviour in various en- 
sembles of stochastic gluon fields developing simple en- 
semble approximation which is grounded on the circum- 
stantial analysis of two-particle correlation function. An 
approximate procedure developed enables to calculate 
the effective functional for the mean energy and to es- 
timate the ground state parameters within the Hartree- 
Fock-Bogolyubov approach. The models with the expo- 
nential and Gaussian behaviours of correlators are ana- 
lyzed in the chiral limit and their parameters are fitted. 
The results obtained are used to estimate the charac- 
teristic region size rg in which the possible processes of 
quark-anti-quark pairing might become significant. This 
size was estimated to herg ^ 1 /40 MeV^^ for the param- 
eters inherent in the NJL model. For the exponential and 
Gaussian correlators this estimate looks like rg ~ 1/150 
MeV-^ 

Besides, we clearly demonstrate the presence of sin- 
gularity in the mean energy functional outside the chi- 
ral limit. Finally, let us emphasize the quark ensemble 
characteristics discussed in the paper are not physically 
observable and in order to make the intelligent conclu- 
sion about the model effectiveness one should explore, 
for example, the meson correlation function. In fact, it 
has been done for the Keldysh model in [sl] and result 
turned out to be quite encouraging. Despite the singular 
character of the mean energy of the system the meson 
observables are finite, quite recognizable and comparable 
with the energy scale coming from an experiment. 
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VI. APPENDIX I 

Here dealing with the simple quantum mechanical ex- 
ample we demonstrate the difference between the descrip- 
tion based on the averaging of density matrix and the 
approach in which the averaging of the generating func- 
tional is used (considering the model system in a real 
time). 

1. Let us suppose the particle described by the station- 
ary Hamiltonian Hq is also affected by the time depen- 
dent force f{t). Then the particle state is circumscribed 
by the Schrodinger equation 

iii , H = Ho + V , V = f X , 
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where ^' is the corresponding wave function and we 
search the perturbative solution as 

\I/ = + + . . . 

We expand the wave function of the zeroth order in the 
eigenfunctions of Hamiltonian Hq as 

= c„ e-*^" , Ho yjn = A„ ^„ , 

the constants c„ are defined by the initial condition here, 
and the next perturbative orders are calculated in the 
following form 

where the coefficient dn^ is determined by the iterations 
as 

dip{t) = -I {^^xi^m) f dr fir) 
Jo 

and (V'„ xtprn) stands here for the matrix element over 
the eigenfunctions of Hq, j = 1,2, ... , dn^ = c„. The 
energy operator after averaging over the final state ^' 
can be presented in the form of trace Tr {Hp} with the 

pseudo-density matrix p = \I>. In general, an energy 
being averaged over such a pseudo-matrix will be time de- 
pendent but at analysing the quasi-stationary processes, 
for example, it might be useful to study its averages in 
time as reads 

Tr {Hp} ^ [ dt Tt{Hp}/T. 
Jo 

For the sake of clarity we suppose for the force mean 
value that / = 0. Then a nontrivial contribution into 
the interaction mean energy comes from the cross terms 
of the zeroth and first orders of perturbation expansion 

^(0) ^(1)^ ^(0) ^(1) and we have 

^ * * 

Tr {Vp} = i CkCm {lpkX1p„){tpnX1p^) X 

(38) 

X e'C^"-^*")* / dr f{T)f{t) e<^''-^"'>'' +C.C . 
Jo 

At estimating the impact of the stochastic force its 
contribution may be factorized with a help of the cor- 
responding correlation function f{T)f{t) ~ F{t — t) 
if the characteristic frequencies A„ are smaller then the 
stochastic ones. In the particular case of the 'white noise' 



(when the profile function F has the (5-function shape) 

the time dependence in the intermediate states ip^^ , ipn in 

Eq. psp disappears (see the corresponding exponentials 

depending on A„). Due to the assumed completeness of 

* 

eigenvalues basis of Hq, i.e. J2n \^n){ipn\ — Eq. ([55)1 
may be presented as 

Tr {^p^^)} ^ jcfe *c,n i^pkx^ i^JP e*(^™-^^)*-f c.c . 
It allows us to conclude that the resulting averaged final 
state density matrix is weighed with the effective 'poten- 
tial' of form /2 x^. The similar results can be received 
in the next perturbative orders. The cluster decomposi- 
tion of stochastic exponential is practical to demon- 
strate that the same results for the effective 'potential' 
of interaction take place for a 'white noise' at averaging 
the generating functional (^) (as it is claimed in the first 
section of this paper). In general consideration there ap- 
pears a certain nonlocal 'potential' and its properties are 
dependent of the system state. 

2. Now let us turn to the description in terms of a den- 
sity matrix only. It is defined by the following equations 

i p = H p- p H' , (39) 

and the density matrix is dependent on the coordinates 
and time p{x, y; t). The operator H is acting on the coor- 
dinate X and the operator H' is acting on the coordinate 
y. In the zeroth order of perturbative expansion we have 

= c„™ e*^"* (y) e-*^-Vm(x) , 

where Cnm is a hermitian matrix which is defined by the 
initial data. In the first order of the perturbation series 
we present the solution in the following form 

= e*^"* V-™ (y) e-*^"V™(x) . 

It is possible to have for the matrix d'^^^ the representa- 
tion as 

d^nl = c„fc (V-™ x^k) f dr fir) e^i^"-^")- + 

Jo 

+i Cfc™ (V-fe y^n) I dT fir) e'i^"-^"^^^ . 
Jo 

Then for the density matrix p^-^"^ the following form ap- 
pears (in order to get it we have to interchange the indices 
m and n in the second term) 
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= - ic„fe e^^^"-^™)* V„ (y) ^m{x) / dr /(r) e'(^"-^'=)" + 



+ ickn [i^k y^m) e'^^"-^")* V™ (y) V'n(x) / dr fir) e^i^"-^-)- 



r 



Due to the hermitian property of density matrix we 

have Ckn —Cnk the second term is complex conju- 
gate with the first one at the coinciding arguments 
X = y. Then calculating the mean interaction energy 
t)\y^x\j we are convinced that the result 
is identical to what we found out at the beginning of this 
section. 

Let us consider the solution for density matrix in the 
operator form (and without specifying the basis func- 
tions) as 

where the matrix p is determined by the solution of fol- 
lowing integral equation 

-pit) = -i I dTe''^"°-">f{T){x-y)e''^"'o~"°^^piT)+p{0). 



and effective interaction Hamiltonian is given by the op- 
erator expression as 

Tt{Hp} = Tr|[iJo + fit) x] i-i) e'^"o-Ho)ty^ 
X f dT e'("o-H^)r j(^) p^^^y.r)\y^^ 



3. Further we analyse some details of the particular 
exercise which admits of receiving the overt expressions 
and consider the forced oscillations defined by the Hamil- 
tonian 



Ho 



1 d^ moj 
2rndx^ ^ 2 



The continual integral is exactly calculated for this ex- 
ercise [igj and the presentation of pseudo-density matrix 
which we are interested in looks like 

'<P iy2,t2)'4>ix2,t2) = // dxidyi if (2/2,^2; J/i,ti) X 



(40) 

X K{x2,t2]xi,ti) V {yi,ti)Lp{xi,ti) , 

where ip{xi,ti) is an initial state. The transformation 
kernel is expressed by the overt formula like 



K{x2,t2;xi,ti) = 



^ 2ni sinwT 
where the action is given as 



1/2 



(41) 



S 



2 sinwT 



[cos LLlT{x2+Xl)-2x2Xi+2x24>2 + '2xi(j)i-F], 



and the phase factor have the following form 
1 Z"*^ 

'^1 ^ ZTT / drfir) sin Lj{t2 - t) , 

02 = / (iT/(r) sin w (r - <i) , 

■mi^ Jti 

T ~ t2 — ti. The term F depends on the time parame- 
ters only and is immaterial because it is cancelled in the 
exponential exponent of pseudo-density matrix. We in- 
troduce the new variable a;2 = ^2 + 0i and transform the 
exponential exponent in the K kernel as 

C0SL0T{x'^ + x\) — 2X2X1 + 2X202 + 2xicf)i ~ 

= cosLoTixl^ + x\) — 2x2Xi + 2{coslljT(I)i + (j)2)x2 + 

+ cos UjT4)1 + 20102 ■ 

The similar transformations should be done in the kernel 

K with the variable 2/2 = ^2 + 0i- The formulae take 
more convenient form if we introduce the auxiliary factor 
01 = (costi;T0i +4)2) /sin ujT which can be written down 
by the simple transformations in the following form 



1 



t2 



dTf{T) C0SUj{t2 — t) 



Making use the well known representation of the kernel 
K for non-perturbated oscillator (/ = 0), see [l9l |. 

Ka{x2M\Xl,ti) = 2^^ '"^"^ 'll>n{x2)i^n{xi). 

It is easy to understand that the important terms of 

pseudo-density matrix kernel KK at 2/2 X2 are rep- 
resented in the similar form 

00 

KK - ^e^^"^ V4y2) k ivi) e~™"^^^^ x 

n=0 



X e-™-i^^ e {X2)i^m{xi) \y^=^^. 



The following matrix element 



dx2 Mh) e-^"'^^^«^ H{X2) e™"^^^^ ^„ (2^2)1,,=., 



will be faced at calculating the mean energy. Now moving 
the exponential e*™'^'^^^^ to the left and changing the 
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variable X2 = X2+4>i in the Hamiltonian H{x2) we obtain 
the representation 

e-'^^'^^^y^ H{x2) e'""^^'^' = Ho{x2) - tu^i i-^ + 

(42) 

+ (mw2(/)i + j)x2 + ^(0? + ^\) + 

which allows us to see that the mean energy calculated 
over the final state is expressed by the diagonal elements 
and matrix elements of coordinate and momentum as well 

/°° * °° * 

dx2 K H {X2) ^ Hn^n 1pn{xi) l/'n (^l) + 

-°° ri=0 
00 

>^ 

+ J2 ^n,n-l e'^^^lPnixi) (yi) + 

00 

* 



where 



n=0 



1 \ mo;^ 
"+2 " + — 



Hn,n-1—H 



n-l,n= [mUJ^{i(t)i + (t)i) + f] (^^) 



1/2 



4. Now we would like to analyse the example of oscil- 
lations initiated by a periodic perturbation defined as 

f{t) = i^sinm . 

Then the phase factors develop the following form 



muj (h-i — F 



mu! (pi = F 



; (ri sin LuT — LU sin ilT) 
r (cosrJT-coswT) , 



(43) 



(for the sake of simplicity we take the parameter as 
ti = 0). As in the limit of classical mechanics these 
expressions include a resonance behaviour at coinciding 
the external frequency Q and oscillator frequency w, and 
in the resonance vicinity (f2 = cj + e with the small devia- 
tion e from the oscillator frequency) the motion behaves 



as the beats, i.e. the small oscillations with the frequency 
uj and large amplitude. Now we are going to resolve the 
corresponding classical equation 

Xc + io^Xc = —f I'm , 
with the initial conditions as Xc(0) = 0, ic(0) — 
F n sin tot ~ LU sin fit 



Comparing Eqs. (|43p and Xc, Xc we are able to express the 
phase factors (pi and (f>i as the classical coordinates (j>i — 
Xc and velocity (pi = Xc/tu- In particular, the correction 
to the diagonal element of effective Hamiltonian (|42|) can 
be presented in the following form 



fPi = ys;^ 



-Xc + fXc 



Averaging the mean energy with pseudo-density ma- 
trix we get the quadratic form as 



(44) 



n,m— 



with the coefficients c„ defined by the initial state and 
normalized as ^ |c„p = 1. In the considered situation 
of periodic force acting for very long (unlimited) time it 
becomes clear the value of mean energy received is time 
dependent (analogously to the classical description) and 
it means a certain asymptotic value for mean energy (as 
for other observables) simply does not exist [isj . Physi- 
cal meaning of this fact appears quite transparent. The 
quantum system is carrying out the repeated transitions 
to the upper levels of excited state and back (these transi- 
tions are controlled by pseudo-density matrix) eventually 
resulting in some quasi-stationary process which can be 
pithily characterized by some observable values averaged 
in time. Thus, the averaged magnitude of diagonal el- 
ement of pseudo-density matrix 7?„^„ — jj" dt iJ„_„/r 
takes the form 



F2 



3cj2 



iF' 



sin2f7r 



m - ijj'^Y 8m r^^- 

— 2aj sin(r2 - 



F' 



2mu;2 ^2 _ ^2 



flT 



n + 2uj sin{Q + uj)T' 
{fl~uj)T n + u! {n + uj)T 



Then it is not difficult to see that asymptotically a quasi- the whole can be set in, indeed, and now the question of 
stationary regime of quantum ensemble oscillations as interest is to determine the minimum of functional (j44p 



16 



which corresponds to some effective ground state of the 
system while under the external influence. The effective 
Hamiltonian (l42]) 



if 



(x) = e-™'^'^i2' H{x) e™'^'^!^ , 



can be presented by using the classical variables Xc, ic 
in the following form 



He f fix) 



iP + Pc) , muj'^ 



(x + Xc) +fix + Xc), 



where pc = mxc- This quantity (at the certain condi- 
tions) may be treated in such a way that it is practical 
to search the ground state with the biased coordinate Xc 
and momentum pc- Certainly, the treatment of excited 
states turns out the nontrivial problem in this situation. 
At every time moment the pseudo-density matrix is a 
pure magnitude because the equality = p is identi- 
cally valid. However, it is possible to estimate the purity 
degree of trail quasi-stationary state = Tr {p^}, with 
the time averaged density matrix and to find such states 
which allow us to develop a description close to one in 
the terms of the Schrodinger equation. 

The density matrix formalism is very practical in more 
general situations, for example, at studying the influence 
of other quantum ensembles on a particle. It is very 
actively discussed and developing (being often quite far 
from our concerns) [34] but our purpose here was to illus- 
trate the difference in describing a system with averaging 
a density matrix and averaging a generating functional. 



VII. APPENDIX II 

The standard way to formulate an effective theory is 
to use the path integral formalism. In order to transit to 
such a description we should construct the corresponding 
Lagrangian action density from the effective Hamiltonian 

(dl 



C = q{iji_,di_,+im)q- G'q t°'-ff_,q J dyli^^{x- y)q't"'-/,,q' , 

(45) 

where q = q(x,i), q = g(x,t), q' = <7(y,<), q' = q(y,t). 
For the highest order in Nc the sum of colour group gen- 
erators looks like J2a=i ^ Kj^ki ~ \ SaSkj- For the sake of 
simplicity we consider the correlator of the following form 
only Ifiu{:s. — y) = <5^,iy /(x — y). Using the Fierz trans- 
formation 7p (g) 7p = 1 (8) 1 + i75 (8) «75 - 57m (8) 7m - 
5 7m 75 'S' 7^75' ^^'^ holding only the scalar contribution 
we receive in the mean field approximation the following 
effective Lagrangian density 



C = q{l^^,d^,^im)q~G' j dy /(x-y) {qq')q' q. (46) 

The brackets in this expression imply the calculation of 
the corresponding averages. The self-consistency condi- 
tion of approximation which may be formulated as the 



following integral equation 
dq 



iM{p) 



(27r) 



G' I{p - q) Tr 



q + im + iM(q) ' 
(47) 

allows us to calculate the quark mass. Integrating over 
the fourth component of momentum 



dq4 



1 



^ 27r ql + q2 + {m + Miq)y 
1 1 



2 [q2-f (m + M(q))2] 



211/2 



we have 
M(p) = 2G'N, 



dq 

(2^ 



/(p-q) 



m + M (q) 



(48) 

With the correlator corresponding to the NJL model we 
obtain the well known gap equation 



M = 2G'Nr 



dq 



m + M 



(27r)3 [q2 + („^ + M)2]i/2 ' 
For the Keldysh model we have 

m + M{p) 



M(p) = 2G'Nc 



[p2 + (m-FM(p))2] 



211/2 



and remember that /(p) = (27r)'^(5(p). Then it is easy 
to understand that the solution can be presented as a 
function of p{M) for convenient handling. 

In the Hartree-Fock-Bogolyubov approach the follow- 
ing sum over the colour matrices is used X^fii ^ ^fj^'jk = 



^ik and then we have for the quark mass 
- 1 



Mg{p) - 2G' 



3 ^(p-q) — sm6'(g) . 



(2^) 



194 1 



(49) 

Comparing this expression to Eg. ([48)1 it becomes clear 
that the four- fermion interaction constant acquires the 
small correction ^ 1/A^c which is rooted in the mean 
field approximation while the higher order terms in 
are held. The patent formula for the Keldysh model can 
be simply received in the chiral limit. In the mean field 
approximation we have 

A/(p)=[(2G'A^,)2_p2]i/2 
and in the Hartree-Fock-Bogolyubov we receive 

1/2 



M{p) 



2G" 



Ni 



1 



p 



see also Ea. ([5S|) . At to 7^ the momentum dependencies 
of masses are quite different. For the Keldysh model in 
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the mean field approximation at zero momentum, for ex- 
ample, we have M(0) = 2G'Nc whereas in the Hartree- 
Fo approach Mg{Q) — Q. At large momenta the mass 
in the mean field approximation behaves as |M(p)| — > 
2G'Nc m/p and in the Hartree-Fock-Bogolyubov approx- 
imation it is the following |M(p)| (^2G' ^jf^ ^ m/p^, 
see also Fig. [3] However, generally, if one takes an orien- 
tation to the analysis of integral characteristics M{p) the 
results are not so different. The similar relations could 
be obtained for the NJL model as well. Apparently it is 
reasonable to notice here that our analysis of the Hamil- 
tonian Eq.(fT2|) (Lagrangians Eqs. (|45l) . (jiG]) ) is also valid 
for the Lorentz-invariant formulation when the r ('time') 
integration is performed for the quark fields as well (for 
the Lagrangians ([IS]) and the integration over 'time' 
is retained and the formfactor becomes a function of four- 
vector I{x — y)). 

= 9 («7m^m + im) q-G / dy I{x-y) {q q') q' q . (50) 



The selfconsistency condition Ea. p7|) acquires the co- 
variant form. In particular, in the Keldysh model in 
four-dimensional formulation when I{p) = (27r)'' S{p) the 
mass gap equation reads as 



M = AGNr 



m + M 



p^ + {m + My 



It allows us to conclude that a quark never comes on the 
mass shell because 



(to + Mf = ANc G 



M 



M 



> 



This feature has already been noticed in Ref. j35i]. The 
similar behaviour has been also observed in the analytic 
models of confinement [36]. Meanwhile, an absence of 
bound states in the four-dimensional Keldysh model (un- 
like the Keldysh model with three-dimensional formfac- 
tors) is its shortage. There appears the additional in- 
tegration over the fourth component of auxiliary four- 
momentum I in Eq. 



(2^ 



2^ 



lih) 



\Pi 



1941 



i h 



where 1(^4) is the respective part of the formfactor. In 
particular, for the four-dimensional Keldysh model with 
/(Z4) = 2-K 6(14) the mean energy functional can be pre- 
sented in the following way 



W{m) 



dp 

(2^ 



\p4 (1 



cos( 



-G 



1 



|P4P 2|p4 



sm ( 



■ cost 



(51) 



The singularity revealed in three-dimensional Keldysh in the four-dimensional consideration, 
model manifests itself as weaker (logarithmic only) one 
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